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ABSTRACT 

In the first week of December 2009, the blazar 3C 454.3 became the brightest high energy 
source in the sky. Its photon flux reached and surpassed the level of 10^^ ph cm~^ s~^ above 
100 MeV. The Swift satellite observed the source several times during the period of high 7- 
ray flux, and we can construct really simultaneous spectral energy distributions (SED) before, 
during, and after the luminosity peak. Our main findings are: i) the optical. X-ray and 7-ray 
fluxes correlate; ii) the 7-ray flux varies quadratically (or even more) with the optical flux; 
iii) a simple one-zone synchrotron inverse Compton model can account for all the consid- 
ered SED; iv) in this framework the 7-ray vs optical flux correlation can be explained if the 
magnetic field is slightly fainter when the overall jet luminosity is stronger; v) the power that 
the jet spent to produce the peak 7-ray luminosity is of the same order, or larger, than the 
accretion disk luminosity. 
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1 INTRODUCTION 

Among the Flat-Spectrum Radio Quasar (FSRQ) class of blazars, 
3C 454.3 (z = 0.859, Jackson & Browne 1991) is one of the 
brightest and most variable sources. In April-May 2005 it under- 
went a dramatic optical outburst, reaching its historical maximum 
with R = 12.0 (Villata et al. 2006). This exceptional event trig- 
gered observations at X-ray energies with the INTEGRAL (Plan et 
al. 2006) and the Swift satellites (Giommi et al. 2006). These de- 
tected the source up to 200 keV and in a brightened high energy 
state. This strong optical and hard X-ray activity was followed by 
a radio outburst with about one year delay (Villata et al. 2007). De- 
spite this extraordinary activity, the lack of a 7-ray satellite at that 
time prevented us from concluding that this activity was due to a 
real increase of the jet power, since the integrated optical luminos- 
ity was of the same order of the 7-ray luminosity observed, years 
before, by the EGRET instrument onboard the Compton Gamma 
Ray Observatory (CGRO) satellite (Nandikotkur et al. 2007). After 
the launch of AGILE, 3C 454.3 was seen to be very bright (Ver- 
cellone et al. 2007) and variable in the 7-ray band, being often the 
brightest 7-ray blazar in the sky (Vercellone et al. 2010), even if the 
simultaneous optical states were much fainter than during the 2005 
flare. Many other active phases followed also in summer 2008 just 
after the launch of the Fermi satellite (Tosti et al. 2008; Abdo et al. 
2009a), and in summer-fall 2009, when the source was detected in 
active state all across the entire electromagnetic spectrum leading 
to many Astronomer's Telegrams (ATels) (Gurwell 2009; Buxton et 
al. 2009; Hill 2009; Striani et al. 2009a; Bonning et al. 2009; Vil- 
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lata et al. 2009). This activity climaxed in December 2009, when 
an extraordinary 7-ray flare occurred on December the 2nd, seen 
by AGILE (Striani et al. 2009b), Fermi/LKT (Escande & Tanaka 
2009), Swift/XKT (Sakamoto et al. 2009), Swift/BAT (Krimm et al. 
2009). 

In the case of the 2007 flare it was possible to model its spec- 
tral energy distribution (SED) in different states with a jet whose 
total power was nearly constant in time (Katarzynski & Ghis- 
ellini 2007), but whose dissipation region was located at differ- 
ent distances -Rdiss from the central power-house (Ghisellini et 
al. 2007). In this model, the magnetic field energy density scales 
as Ub oc (r_Rdiss)~'^, where F is the bulk Lorentz factor. Con- 
versely, the radiation energy density (7cxt produced by the broad 
line clouds scales as and is independent of -Rdiss as long as the 
dissipation region is within the broad line region (BLR). Even if 
the latter is most likely stratified, the contribution of the Ly-a line 
to the BLR luminosity is dominant, and we identify the size of the 
BLR with the shell mostly emitting the Ly-a photons. As a conse- 
quence, Us/Uext ~ Lopt/Lj oc i^diss'^"*- Even with a constant 
and dominating and a constant F, we can have large variations 
of the optical flux by changing -Rdiss - These ideas were borne out in 
Ghisellini et al. (2007), immediately after the AGILE detection of a 
very bright 7-ray state accompanied by a moderate optical flux: the 
optical synchrotron flux can vary by a large amount even if the jet 
continues to carry the same power and always dissipates the same 
fraction of it. This assumes and implies, however, that the bulk of 
the electromagnetic output is in the high energy (inverse Compton) 
component, that in this case should remain quasi-steady. 

Instead the large outburst seen in December 2009 witnessed 
a dramatic increase of this component. In this case it is very likely 
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that the power dissipated by the jet indeed increased, and it is hkely 
that this is the consequence of an increase of total power carried by 
the jet. The luminosity reached in 7-rays exceeded L-y ~ 3 x 10^^ 
erg s~^, a factor ~30 larger than the EGRET luminosity, and also 
larger than the 2007 AGILE state. The goal of our present study is 
first to characterise the emission properties of 3C 454.3 during the 
November-December 2009 period, then to find the physical prop- 
erties of the jet during the flare, and finally to compare those with 
the luminosity of the accretion disk. In fact, as we shall see, the 
power of the jet of 3C 454.3 in the first week of December 2009 
likely surpassed the luminosity of the accretion disk. 

We use Hq = 70 km s"^ Mpc"^ and VLm = 0.3, J^a = 0.7 
and use the notation Q — 10^ Qa; in cgs units. 



2 DATA & ANALYSIS 

2.1 Data sets 

This study is focused on the analysis of the multi-wavelength 
(MWL) data gathered by the Fermi and Swift satellites in the period 
between November 1 and December 14, 2009. From these data we 
obtained day-scale resolved light curves along the whole time span, 
in the 7-ray, X-ray and optical-UV band. Moreover, for some se- 
lected days we derived the spectral information needed to study the 
simultaneous MWL SED. These data sets are labelled according to 
the observation date. 

In order to fully characterise the evolution of the SED across 
the whole dynamic range of 7-ray flux observed by Fermi, we also 
built a "low" SED, deriving a high-energy (HE) spectrum of 3C 
454.3 from an earlier period of observation when the source was 
significantly weaker in 7-rays, and matching it with Swift data from 
the same period. 

Finally, in order to better investigate the UV portion of the 
SED, we matched an archival UV spectrum of 3C454.3 observed 
by the GALEX satellite with contemporary optical-UV data taken 
by Swift. 

2.2 Fermi/LAT 

We retrieved from the NASA databas^B the publicly available data 
from the LAT 7-ray telescope onboard of the Fermi satellite (At- 
wood et al. 2009). We selected the good quality ("DIFFUSE" class) 
events observed within 10° from the source position (J2000) (a, 
(5)=(22h53m57.7s, +16°08'53.56"), but excluding those observed 
with zenith distance of the arrival direction greater than 105°, in 
order to avoid contamination from the Earth albedo. We performed 
the analysis by means of the standard Science Tools, v. 9.15.2, in- 
cluding Galactic and isotropic extragalactic backgrounds and the 
P6 V3 DIFFUSE instrumental response function. For each data set 
(and, in case of spectra, for each energy bin) we calculated the live- 
time, exposure map and diffuse response of the instrument. Then 
we applied to the data an unbinned likelihood algorithm (gtlike), 
modelling the source spectrum with a power law model, with the 
integral flux in the considered energy band and photon index left as 
free parameters. The minimum statistical confidence level accepted 
for each time (or time and energy, in the case of spectra) bin is TS = 
10 where TS is the test statistic described in Mattox et al. 1996 (see 
also Abdo et al. 2009b). Quoted errors are statistical, at I cr level. 
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Systematic flux uncertainties are within 10% at 100 MeV, 5% at 
500 MeV and 20% at 10 GeV (Rando 2009). 



2.2.1 Light curve 

In the upper panel of Fig. [T] we report the integral 7-ray flux above 
100 MeV, calculated in 1-day time bins for the period from Novem- 
ber I to December 14, 2009. The flux level is always above 10^® 
cm^^ s^^, the clear signature of a coherent 7-ray active phase. On 
December 2 the flux reached the extreme value of (21. 8± 1.2) x 
I0~^ cm~^ s~^, the strongest non-GRB flux measured by a 7-ray 
satellite. The corresponding observed 7-ray luminosity is of the 
same order of PKS 1502+106 in outburst (Abdo et al. 2010a). 



2.2.2 Spectra 

The exceptional flux level of 3C 454.3 in late 2009 allowed for 
the first time to derive day-scale 7-ray spectra with good quality: 
usually an inconceivable task for space-based instruments due both 
to their small collection area and the weak flux of the high-energy 
sources. This opens new opportunities for the study of the MWL 
SED of this FSRQ on short time scales, hopefully allowing deeper 
insight into the engine and the emission mechanism. 

In this study we use LAT spectra produced with data from 
only a few single days, selected within the period covered also by 
Swift observations. We have chosen to construct the SEDs in order 
to cover the whole flux range of the campaign. We therefore se- 
lected November 6, showing the lowest UVOT fluxes of the whole 
period; November 27, a day in the rising phase of the flare, with 
intermediate 7-ray flux, and the three days (December 1, 2 and 
3) around the peak, in order to better investigate the evolution of 
the SED around this extreme phase. For the 7-ray spectrum of 
November 6 we decided to add the data from the following day 
(November 7), as it showed compatible 7-ray flux level, to im- 
prove the signal to noise ratio (S/N) of the spectrum. Finally, we 
selected from the 18-months daily light curve (see e.g. Tavecchio 
et al. 2010) a longer time-span (15 Ms), starting from 3 December 
2008 to 26 May 2009, when the source was at the lowest flux levels 
{F>ioo Mcv ~ 5 X 10~^ cm~^ s~^) since the beginning of LAT ob- 
servations. The derived HE spectrum allowed to explore the whole 
available dynamic range of 3C 454.3 7-ray emission. For each data 
set, we split the LAT observed band in logarithmically spaced en- 
ergy bins, and applied in each the unbinned likelihood algorithm. In 
addition to the request that TS > 10, we further requested that the 
model predicted at least 5 photons from the source; when needed, 
we merge in a single, wider bin the high energy bins that did not 
fulfil this last request. These quality criteria are similar (and some- 
how more restrictive) to those adopted by Abdo et al. (2010b). For 
each spectrum, we crosschecked this analysis performed in bins of 
energy with an unbinned analysis on the whole O.I-IOO GeV band, 
where we modelled the source with a broken power law model, in 
agreement with Abdo et al. (2010b). The results of the two methods 
were consistent within errors, even if the analysis in bins of energy, 
that does not exploit the full information of the event distribution 
across the whole LAT bandpass, leads to wider error bars. We re- 
port in Table[T]the parameters of the model, while the SED points 
are plotted in Fig. [4] 
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2.3 Swift 

We analysed the data from observations performed between 
November 1 and December 14, 2009. We also analysed data from 
late December 2008, contemporary to the lowest 7-ray state of 3C 
454.3 seen by Fermi. Finally, UVOT data from a pointing in Octo- 
ber 2008 were analysed, as we wanted to compare this observation 
with a roughly simultaneous GALEX archival UV spectrum of the 
source (see i]3.2n . 

The data from the two narrow field instruments onboard Swift 
(Gehrels et al. 2004) have been processed and analysed with 
HEASof t V . 6.8 and with the CALDB updated on December 
30, 2009. 



2.3.1 X—Ray Telescope 

Data from the X-ray Telescope (XRT) (0.2-10 keV, Burrows et 
al. 2005) were analysed using the xrtpipeline task, set for 
the photon counting or window timing modes and having selected 
single pixel events (grade 0). In the case of the light curve, we 
converted the observed counts rates to fluxes assuming a common 
r = —1.6 photon index and with absorption fixed at the Galactic 
value N^""^ = 6.6 x 10^" cm~^ derived from Ab = 0.46 (Raiteri 
et al. 2008). The resulting light curve is shown in the middle panel 
ofFig.E 

We extracted the spectra in order to build SED simultaneous 
to FermifLAT by using default regions. Data were rebinned in order 
to have at least 30 counts per energy bin. Power law models have 
been fitted to the spectra, except for two cases, the "low" state and 
the spectrum from December 2, when a broken power law model 
fitted the data significantly better ( ftest > 99%). In Tab.|2]we report 
the spectral parameters for the observations used for modelling the 
SED. 



2.3.2 Ultraviolet— Optical Telescope 

Data from the UV telescope UVOT (Roming et al. 2005) were anal- 
ysed by means of the uvotimsum and uvotsource tasks with 
a source region of 5", while the background was extracted from a 
source-free circular region with radius equal to 50" (it was not pos- 
sible to use an annular region, because of a nearby source). A total 
of 19 UVOT observations have been performed and analyzed be- 
tween November 1 and December 14. The extracted fluxes are plot- 
ted in the lower panel of Fig. [T] The Galactic extinction has been 
corrected assuming Ab = 0.46 and Rv ~ Av/E{B — V) = 3.1 
(according to Raiteri et al. 2008), and the extinction law of Cardelli 
et al. (1989). The observed magnitudes from the 5 pointings used 
for the MWL SED are reported in Table [3j uncertainties include 
systematics. The observation made on October 1, 2008, roughly 
contemporary to the GALEX archival spectrum (see i]2.4t is also 
reported in the table. The UVOT data set used for the modelling 
of the "low" state is built with pointings made between December 
25, 2008 and January 1, 2009. In this case the reported magnitudes 
are the mean values over the observations, while the uncertainties 
are calculated as cr = V^tat+^'ftd' where o"stat is the mean sta- 
tistical uncertainty on the single observation and (Tstd is the sample 
standard deviation of the magnitude distribution, and keeps into ac- 
count the variability of the source. 
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Figure 1. Light curve above 100 MeV (top), in the [0.2-10 keV] X-ray 
band (middle panel) and in the different UVOT filters (bottom panel, as 
labelled). Note the logarithmic y axes, and their different scales. Yellow 
lines mark the days for which simultaneous multi-wavelenght SED were 
derived and modelled. 



2.3.3 Burst Alert Telescope 

Swift monitors the sky in the hard X-ray band by means of the 
Burst Alert Telescope (BAT), a coded mask telescope sensitive in 
the 15-150 keV band (Barthelmy et al. 2005). We retrieved from 
the HEASARC hard X-ray transient monitor databas^ the daily 
sampled light curve of 3C 454.3 in the 15-50 keV band, and 
converted the count rates to fluxes assuming as reference the BAT 
count rates and flux of the Crab Nebula in the same band. The cor- 
responding vFu values are plotted in Fig.|4] 



2.4 GALEX 

The Galaxy Evolution Explorer (Martin et al. 2005) is a NASA 
satellite, in flight since 28 April 2003 and performing an all-sky 
survey in the far UV (FUV, ~ 154 nm) and near UV (NUV, ~ 232 
nm) with resolution R = A/ A A ~ 200 in NUV slit-less spec- 
troscopy (Morissey et al. 2007). We retrieved from the MASlff] 
database a publicly available spectrum of 3C 454.3, obtained from 
3 exposures taken between September 30 and October 2, 2008. This 
spectrum is plotted in Fig. [3] and its relevance is discussed in Wi.2\ 



^ http://heasarc.gsfc.nasa.gov/docs/swift/results/transients/ 

Multimission Archive at the Space Telescope Science Institute, 

http : / /galex . St sci . edu/GR4 / 
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OBS date 






r2 


77 

-'-'break 


\i . L — - lUU "-jC V 


TS 


dd/mm/yy 


ks 






GeV 


10-6^ 




low 


5.9x10^ 


2.33±0.03 


3.37±0.27 


2.00±0.38 


0.61 ±0.02 


4687 


06/11/09 


61 


2.66±0.21 


3.80±1.20 


1.10±0.57 


2.10±0.30 


218 






2.78±0.16 






2.15±0.30 


218 


27/1 1/09 


37 


2.13±0.10 


3.17±0.54 


1.94±0.56 


7.58±0.68 


900 


01/12/09 


44 


2.33±0.11 


3.15±0.38 


0.97±0.34 


10.76±0.87 


1100 


02/12/09 


39 


2.37±0.06 


3.17±0.46 


2.70±0.26 


21.6±1.2 


2497 


03/12/09 


40 


2.28±0.06 


3.42±0.87 


5.86±0.36 


15.94±0.98 


1915 



Table 1. Results of the 7-ray analysis on Fermf'-LAT data; only the data set used for the MWL SED are reported. For each data set, the effective source 
on-time ton, the photon indices Fi and ^2, the break energy Bbrcak, the integral flux in the 0.1 - 100 GeV band Fg i„if)o GcV, and the TS value of the 
unbinned likelihood analysis are reported. Uncertainties are statistical only, at 1 cr level; systematics on flux measurement are within 10 % at 100 MeV, within 
5% at 500 MeV and within 20% at 10 GeV. For the 06/1 1/009 dataset, the relative weakness of the flux makes perfectly adequate a simple power law model 
which is also reported for reference. 



OBS date 


obslD 


^cxp 






-^brcak 


/o 


-'^o.a-iokcv 


Lo.2-10 kcV 


(dof) 


dd/mm/yy 




ks 






keV 


10-2 .T\ „ 

cm^ s kc V 


jO_ll_org 

cm^ s 


lO*^ erg/s 




low 

06/11/09 
27/11/09 
01/12/09 
02/12/09 
03/12/09 


see caption 
00035050070 
00035050073 
00035050074 
00035050075 
00035050076 


9.0 
1.0 
0.9 
1.1 
1.2 
1.0 


, ,0+0.30 
^■^-'-0.24 

1.52±0.07 
1.42±0.06 
1.43±0.04 
0.98l°- 
1.43±0.04 


, ,0+0.30 

'■"-0.24 

1.56±0.07 


'■^"-0.36 
1.2 ±0.2 


0.15±0.01 
2.2± 0.2 
3.0 ±0.2 
4.5 ±0.3 
2.2 ±0.2 
5.2 ±0.3 


I. 1 

6.6 

II. 
16.0 
17.0 
18.0 


0.34 
2.1 
3.1 
4.5 
4.8 
5.2 


0.89 (43) 
0.96 (38) 
1.13 (53) 
1.08 (87) 
0.98 (103) 
0.88 (89) 



Table 2. Synopsis of the results of the X-ray analysis performed on SwiftfXKY data. Only the data sets used for modelling the MWL SED of 3C 454.3 are 
reported. For each one of these the exposure time tcxp, the photon indices Fi and F2, the break energy -Ebrcak, the normalisation factor of the differential 
spectrum at 1 keV /q, the observed (not de-absorbed) flux in the 0.2-10 keV band Fo.2— 10 koV, the intrinsic luminosity in the same band Lo.2-10 kcV and 
the goodness of fit statistic (reduced x^) with the coiresponding number of degrees of freedom (dof) are reported. The "/ow" data set results from the 
addition of observations made between 25 Dec 2008 and 01 Jan 2009 (obslD from 00090023002 up to 00090023008) while the other data sets come frorn 
single pointings and are labelled with the date of observation. Spectra have been de-reddened according to the = 6.6 X lO^" cm— ^ value derived from 
Raiteri et al. (2008). The errors on the observed fluxes are within the 10% level (at 90% CL). 



3 OPTICAL AND X-RAY VS 7-RAY FLUX 

Fig. [2] shows the 7-ray flux as a function of the flux in the UVOT 
uvw2 and b filters and as a function of the XRT X-ray flux. From 
the uvw2 flux we have subtracted the contribution of the accretion 
disk (0.55 mjy) and the tail of the Ly-a line (contributing for 0.094 
mjy) intercepted by this filter (see Fig.O. For the b fluxes, we have 
subtracted a disk contribution of 0.64 mJy. 

The total time interval is the same shown in Fig.[T] i.e. about 
1.5 months. The 7-ray flux correlates with both the UV-optical 
and X-ray fluxes, and the correlation is more than linear: for a total 
variation of a factor ~20 of Ky, the optical-UV and the X-ray vari- 
ability amplitude is a factor ~3^. The vertical grey stripe indicates 
the flux contribution in the m>w2 UVOT filter due to the accretion 
disk emission according to our model plus the contribution from 
the tail of the Ly-a line intercepted by the uvw2 filter. 

Formally, a least square fit assuming a relation log = 
mlogFuvw2 + q yields m = 1.567 and q = —5.026 (these are 
the values for the bisector of the correlations y = mix + q\ and 
X = may + 52), with a chance probability P = 6.4 x 10"^. 

The correlation with the X-ray flux, as illustrated in Fig. 2, 
seems to suggest the presence of an X-ray constant component 
contributing at a level of ~ 5 x 10~^^ erg cm-^ . On the 
other hand, the plotted data refer only to the Nov. - Dec. 2009 



campaign, while the "low" state (shown in Fig.|4ll shows that there 
are states with lower X-ray and 7-ray fluxes (see Tab. [T] and Tab 
[2l(. The suspected candidate for this behaviour is the SSC com- 
ponent during the Nov.-Dec. 2009 period, as we will discuss in 
§5.1. Applying a least square fit to the data of Fig. [2] we find 
log = 2.267 log Fx + 17.512 (again these are the bisector val- 
ues) and a chance probability of P = 5 x 10~^. 

The correlated variability in the three bands strongly suggests 
that the corresponding fluxes originate in the same region and by 
the same population of electrons. This supports the "one-zone" 
models. 

The fact that F-y varies more than linearly with the -optical- 
UV and the X-ray fluxes is interesting and, at first sight, surpris- 
ing. Indeed the optical-UV is likely synchrotron emission, while 
the X-ray flux, having a shape much harder than the optical spec- 
trum, belongs to the same high energy hump as the 7-ray emission. 
We then expect that F-y and Fx vary together, with F-y oc Fx ■ Fur- 
thermore, if the 7-ray emission is produced by inverse Compton 
scattering with broad line photons, produced externally to the jet 
- i.e. it is External Compton (EC) radiation - it is proportional to 
the number of the emitting electrons, as the synchrotron flux. Vari- 
ations of the electron number should then result in Ky oc Fopt . Fi- 
nally, if the synchrotron Self-Compton model (SSC) process con- 



The ^-ray brightest days of the blazar 3C 454.3 5 



OBS date 


obsID 


11 


b 


u 


uvwl 


uvm2 


uvw2 


01/10/08 (GALEX) 


00031216062 


15.42±0.04 


15.92±0.03 


15.24±0.04 


15.41±0.04 


15.50±0.04 


I5.7I±0.04 


low 


see caption 


16.23±0.11 


16.72±0.07 


15.93±0.08 


15.96±0.09 


15.92±0.08 


16.18±0.05 


06/11/09 


00035050070 


15.88±0.06 


16.33±0.05 


15.60±0.05 


15.68±0.05 


15.68±0.06 


15.86±0.05 


27/11/09 


00035050073 


14.77±0.04 


15.25±0.04 


14.58±0.04 


14.85±0.05 


14.95±0.05 


15.13±0.04 


01/12/09 


00035050074 


14.67±0.04 


15.20±0.03 


14.57±0.04 


14.79±0.04 


14.93±0.05 


15.16±0.04 


02/12/09 


00035050075 


14.49±0.04 


15.02±0.03 


14.25±0.04 


14.56±0.04 


14.67±0.05 


14.89±0.04 


03/12/09 


00035050076 








14.58±0.04 


14.75±0.04 


14.97±0.04 



Table 3. Summary of Swift/WOT observed magnitudes. Only the data sets that we used for modelling the MWL SED of 3C 454.3 are reported here, out of 
the total of 19. In addition, we report the UVOT observation that we compared to the GALEX archival spectrum in Fig. [3] The data set used for the "low" state 
includes 7 UVOT pointings, with obslD starting from 00090023002 up to 00090023008. Magnitudes are not de-reddened; errors include systematics. 




[mJy] 

Figure 2. CoiTelation between the non-thermal UV flux (UVOT filter 
uvwl, subtracting the disk and the Ly-a contributions), non-thermal op- 
tical flux (UVOT filter h, subtracting the disk contribution), non-thermal 
X-rays and the 7-ray flux. The grey vertical stripe indicates the contribu- 
tion of the accretion disk plus the tail of the Ly-o line entering in the uvw2 
frequency range. The dashed lines are the linear fit to the 7-ray flux with 
the uvwl and X-ray fluxes (see text). 
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Figure 3. A portion of the UV SED of 3C 454.3 as seen by GALEX in Oc- 
tober 2008 (in black). The red SED points are derived from a contemporary 
UVOT observation. Atz = 0.859 the Ly-a line (Aq = 1216 A) is shifted 
close to the effective frequency of the uvm2 filter, while the neighbouring 
uvwl and uvwl are less sensitive to it. The observer frame frequency of 
the Ly-a line is shown by the dotted vertical line. The effective areas of 
the UVOT filters are also plotted for reference, with thick dotted lines; the 
reference scale is reported on the right vertical axis. 



tributes substantially in the X-ray band, Fx is proportional to the 
square of the number density of the emitting electrons, and then we 
expect Fx oc F^, just the opposite of what observed. A very sim- 
ilar behaviour has been observed for PKS 1502-1-106 (Abdo et al. 
2010a). 

More than linear and even more than quadratic variations of 
the high energy 7-ray flux with respect to the X-ray flux have been 
already observed in TeV BL Lacs for relatively short (i.e. weeks) 
periods of time, as in Mkn 421 (Fossati et al. 2008) and PKS 2155- 
304 (Aharonian et al. 2009), but in these cases the X-ray emission 
belongs to the synchrotron part of the spectrum, while the 7-ray 
flux is probably SSC emission. Even so, it is not easy to explain a 
quadratic relation during the decaying phases of the light curves, 
because this implies short cooling timescales (Katarzynski et al. 
2005), leading Katarzynski & Walczewska (2010) to propose that 
the X-ray flux of the varying component must be diluted by the 
flux produced by another region. 



In conclusion, we have that the 7-ray flux varies more than the 
flux at lower energies both in 3C 454.3 (and PKS 1502-1-106) and 
in the less powerful BL Lac objects, despite the large differences in 
the jet power and in the jet environment leading to different emis- 
sion processes for the 7~ray photons (i.e. SSC for BL Lacs and 
EC for 3C 454.3 and other powerful FSRQs). While these similar- 
ities deserve further studies, we note here that for powerful blazars 
we have very short cooling times of the emitting electrons, even at 
low energies: this at least helps to explain the decaying phases of 
the light curves, that are instead a problem for TeV BL Lacs. For 
the specific case of 3C 454.3, we will suggest in the following the 
likely cause of F.y varying more than Fx and Fopt, in terms of an 
inverse correlation between the dissipated power and the magnetic 
field. 
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3.1 Variability and cooling timescales 

Tavecchio et al. (2010) analysed in detail tlie entire light curve of 
3C 454.3 in 7-rays, including the period of exceptional activity 
studied here. It is found that the minimum doubling timescale is 
between 3 and 6 hours (in our observer frame). This result has been 
also confirmed by Foschini et al. (2010) who studied another flare 
occurred in early April 2010, comparable in flux but observed by 
Fermi-LAT in pointed mode, therefore with denser time coverage 
and enhanced statistics. From Fig. [2] we can see that also in X- 
rays (note the decrease at Dec. 6 and at Dec. 12) the variability 
timescale is ~1 day or shorter. This behaviour has two important 
consequences: 

• The size of the emitting region must be compact, R < 
cU^,5/{l+z) ^ 7xl0^5tvar,6h(5/20) cm. As remarked in Tavec- 
chio et al. (2010), this challenges models in which the dissipation 
region is at tens of parsecs from the black hole (see e.g. Jorstad 
et al. 2010), and instead favours models in which the dissipation 
occurs much closer. 

• The cooling time fcooi ^ tvar for electrons emitting at 7-ray 
and X-ray energies. This favours models in which the dissipation 
occurs in a region where the high density of photons can ensure a 
rapid cooling. Again, dissipation close to the black hole is favoured, 
where the broad line region (BLR) can provide a radiation energy 
density large enough to have fcooi < ^var even for electrons with 
random Lorentz factors of a few. 



3.2 Contribution of the Ly-Q and of the accretion disk 

At a redshift of 2; = 0.859, the Ly-a emission line from the BLR 
is observed at 2260 A, within the band of the uvm2 filter of UVOT 
(effective wavelenght Ae// = 2231 A, Poole et al. 2008). 

In this source the line is bright (rest frame equivalent width 
EW = 74 A according to the observations by Wills et al. 1995), and 
we wanted to check if it could explain the small excess in the uvml 
and uvw2 filters that is visible in the SED (Fig.|4]l. In Fig.|3]we show 
the 1900-2750 A portion of the UV spectrum of 3C 454.3 observed 
by the GALEX satellite in October 2008, together with the SED 
points observed in the UVOT filters and the corresponding trans- 
mission curves. All data have been de-reddened with Ab ~ 0.46 
and Rv ~ 3.1, following Raiteri et al. (2008) and assuming the 
extinction law from Cardelli et al. (1989). A small mismatch is ap- 
parent, but many factors can account for it, such as systematics in 
the flux measurement between the two instruments and the loose si- 
multaneity of the observations performed with the two telescopes. 
Moreover, a quantitative study on this issue is far beyond the scope 
of this paper. Fig.[3]shows that the small bump in the portion of the 
SED observed by the UVOT filters can be partially due to the Ly-a 
emission line shining the uvml, and, at a lower level, the neighbour- 
ing filters, especially uvw2. This excess was well visible in all our 
UVOT pointings despite the larger level of the non-thermal contri- 
bution. Further, we tried to build a simple model for this portion of 
the spectrum, as the sum of three components: 

• the Ly-Q emission line, modelled as a gaussian profile, with 
a and normalisation obtained from a fit to the GALEX spectrum; 

• the blue bump due to a "standard" accretion disk, assumed 
constant (see SQ; 

• a variable power law contribution, representing the syn- 
chrotron emission from the jet. 

With these three components we managed to reproduce the SED 
profile for each UVOT observation we considered. The excess in 



the uvml filter that we attributed to the Ly-a line was compatible, 
within the uncertainties, with a constant contribution. This supports 
Raiteri et al. (2007) when they explain the smaller variability range 
in the UV band w.r.t. the optical and infrared (IR) as due to the 
presence of a blue, stable contribution, becoming more important 
(relative to the steep synchrotron non-thermal continuum) at higher 
frequencies. What we emphasise is the presence of two constant 
components: the accretion disk and the Ly-a line. 



3.3 The SED of 3C 454.3 at different epochs 

In Fig.|4]we show the overall optical to 7-ray SED of 3C 454.3 
at 5 different epochs during November and December 2009 and 
an additional one ("/ow") representative of a low 7-ray state (see 
[J2J. We can see the large variability amplitude of the 7-ray flux, 
whose spectrum is instead remarkably less variable, as noted also 
in Abdo et al. (2010b). Although it is somewhat harder in the "/ow" 
state and on 27 Nov. 2009, it is always characterised by a slope 

> 1 [F{y) oc i^^"], so that the peak of the 7-ray spectrum 
is at energies smaller than 100 MeV On the other hand the X-ray 
0.2-10 keV slope and the Swift IHK^ data points suggest that the 
peak cannot be much below 100 MeV. 

As discussed above, we can see that the optical-UV flux 
varies, but less than the 7-ray flux. At low states, the accretion 
disk becomes "visible" by flattening the spectrum. The presence of 
the accretion disk emission has been noted by Raiteri et al. (2007) 
discussing observations in 2006-2007, that indicate that the disk 
produced a flux comparable with what reported in Fig.|4] 

As remarked before, also the X-ray flux varied much less than 
the 7-ray one during Nov.-Dec. 2009. On the other hand, in the 
"/ow" state, the X-ray flux is much lower (also with respect to the 
optical), and still very hard. 

3.4 The mass of the black hole of 3C 454.3 

Gu, Cao & Jiang (2001) estimated M = 4 x l(f Mq for the mass 
of the black hole of 3C 454.3. They used the size of the BLR JJblr 
and the FWHM of the H/3 fine, equal to 2800 km s~^ (Netzer et 
al. 1995). To derive the velocity of the clouds, they multiplied the 
FWHM value by a factor / = 1.5 (V = /Vfwhm), according to 
McLure & Dunlop (2001). The McLure & Dunlop / = 1.5 factor 
was empirically estimated assuming that the clouds belong to two 
components: one moving randomly and isotropically (that would 
give / = \/3/2) and one distributed in a disk like geometry. To 
derive the black hole mass through M = _RblrV'^G~^, -Rblr 
was found using the relation of Kaspi et al. (2000): 



-Rblr ~ 4.2 x lO' 



H 0.700±0.033 



10^5 ergs- 



(1) 



where (•^^A)gj^oQ^ the luminosity of the continuum at 5100 A 
(rest frame, corresponding to an observed wavelength of 9480 A). 
For high luminosity quasars, the more recent results of Kaspi et al. 
(2007), obtained using the CIV A 1549 line, yield: 

-, 0.55±0.04 



-Rblr ~ 7.8 x 10' 



lO^s ergs- 



(2) 



The observed luminosity at 2510 A (corresponding to 1350 A, rest 
frame) in a period of low state, not much contaminated by the non- 
thermal continuum, was \Lx ~ 4 x 10*® erg s~^ (from the spec- 
trum taken by the GALEX satellite reported in Fig.[3]l. This yields 
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-Rblr ~ 6 X 10^^ cm. The FWHM of the CIV line, according to 
Wills et al. (1995) is 3145 km s^^. Assuming an isotropic velocity 
field, / = \/3/2, the mass calculated through M = Rbi^kV^G'^ 
is M = 3.4 X W^Mq. Assuming / = 1.5, instead, we find 
M = 1.2 X lO^'Af i), in any case smaller than estimated by Gu 
et al. (2001). Since the continuum could in fact be contaminated 
by the non-thermal beamed contribution, these values should be 
considered as upper bounds. Furthermore, given the uncertainties 
in these relations, this estimate can give only a rough indication of 
the black hole mass. In the following, when modelling the SED, we 
will use M = 5 X 10* M©. In the general framework of our model, 
this is the black hole mass ensuring an observed minimum variabil- 
ity timescale of ~6 hours, if the jet opening angle in its conical part 
is 0. 1 rad, as assumed for all other blazars analysed in Ghisellini et 
al. (2010). 



4 MODELLING THE SED 

To model the SED we have used a relatively simple leptonic, one- 
zone synchrotron and inverse Compton model. This model, fully 
discussed in Ghisellini & Tavecchio (2009), has the following main 
characteristics. 

We assume that in a spherical region of radius R, located at 
a distance -Rdiss from the central black hole, relativistic electrons 
are injected at a rate Q{'y) [cm^"^ s^^] for a finite time equal to 
the light crossing time R/c. For the shape of Qij) we adopt a 
smoothly broken power law, with a break at 7b: 



Q(7) = Qo 



1 + (7/7b)-"i+''2 



(3) 



The emitting region is moving with a velocity /3c correspond- 
ing to a bulk Lorentz factor P. We observe the source at the viewing 
angle 9^ and the Doppler factor is (5 = l/[r(l — Scos^v)]. The 
magnetic field B is tangled and uniform throughout the emitting 
region. We take into account several sources of radiation externally 
to the jet: i) the broad line photons, assumed to re-emit 10% of 
the accretion luminosity from a shell-like distribution of clouds lo- 
cated at a distance Eblr = IQ^'' L]/^^ cm; ii) the IR emission 
from a dusty torus, located at a distance _Rir = 2.5 x IQ^^L^J^^ 
cm; iii) the direct emission from the accretion disk, including its 
X-ray corona. Furthermore we take into account the starlight con- 
tribution from the inner region of the host galaxy and the cosmic 
background radiation, but these photon sources are unimportant in 
our case. We instead neglect, for simplicity, the contribution of a 
possible population of intra-cloud free electrons. This contribution 
would be important if the optical depth t of these electrons be- 
came greater than the covering factor of the emitting line clouds, 
i.e. r > 0.1, a rather large value. These free electrons would scat- 
ter and re-isotropize the entire disk flux, providing IR radiation, 
besides optical and UV. The final effect would be to broaden the 
high energy hump, especially below its peak. 

All these contributions are evaluated in the blob comoving 
frame, where we calculate the corresponding inverse Compton ra- 
diation from all these contributions, and then transform into the ob- 
server frame. For simplicity, the numerical code assumes an abrupt 
cut-off for the scattering cross section, equal to the Thomson one 
for ^jhvi^y-a/mcf? < 1 and zero above. This approximation pro- 
vides acceptable results in the calculation of the IC scattering of the 
BLR (e.g. Tavecchio & GhiselUni 2008) and starts to fail only at 
the largest energies probed by LAT, at which scatterings in the full 



Klein-Nishina regime start to be important (see also Georganopou- 
los Kirk & Mastichiadis 2001 ; Moderski et al. 2005 for a discussion 
on the importance of the full Klein-Nishina cross section). 

We calculate the energy distribution A'^(7) [cm~^] of the emit- 
ting particles at the particular time R/c, when the injection process 
ends. Our numerical code solves the continuity equation which in- 
cludes injection, radiative cooling and pair production and re- 
processing. Our is not a time dependent code: we give a "snapshot" 
of the predicted SED at the time R/c, when the particle distribution 
A'^(7) and consequently the produced flux are at their maximum. 

For 3C 454.3 the radiative cooling time of the particles is 
short, shorter than R/c even for low energetic particles. This im- 
plies that, at lower energies, the A'^(7) distribution is proportional 
to 7 while, above 7b, N('y) oc 7^(^2+1) jjjg electrons emit- 
ting most of the observed radiation have energies 7peak which is 
close to 7b (but these two energies are not exactly equal, due to the 
curved injected spectrum). 

The accretion disk component is calculated assuming a stan- 
dard optically thick geometrically thin Shakura & Sunjaev (1973) 
disk. The emission is locally a black body, with a temperature pro- 
file given e.g. in Frank, King & Raine (2002). Since the optical-UV 
is the sum of the accretion disk and the jet non-thermal compo- 
nents, there is some degeneracy when deriving the black hole mass 
and the accretion rate, even in the lowest state, for which the disk 
is best visible. On the other hand we can take as a guideline the 
value of the black hole mass derived in §4.1, implying then a very 
large accretion rate. We also note that in modeling the optical-UV 
region one should also consider other thermal components, such as 
the other emission lines and the Balmer continuum from the BLR. 
However, all these components are expected to be much less bright 
than the Ly— a emission Une, included in the fit. It is also possible 
that the disk emission spectrum is different than that expected by 
the Shakura & Sunjaev model: however we use the latter, built on 
robust theoretical basis. 

We model at the same time the thermal disk (and IR torus) ra- 
diation and the non-thermal jet-emission. The link between these 
two components is given by the amount of radiation energy den- 
sity (as seen in the comoving frame of the emitting blob) coming 
directly from the accretion disk or reprocessed by the BLR and the 
IR torus. This radiation energy density depends mainly on J?diss, 
but not on the adopted accretion rate or black hole mass (they are 
in any case chosen to reproduce the observed thermal disk lumi- 
nosity). 

By estimating the physical parameters of the source we can 
calculate the power that the jet carries in the form of radiation (Pr), 
magnetic field (Pb), relativistic electrons (Pc) and cold protons 
(Pp) assuming one proton per electron. These powers are calcu- 
lated according to: 



Pi = nR^r'^cUl 



(4) 



where U' is the energy density of the ith component in the comov- 
ing frame. 



4.1 Constraints on the model parameters 

We briefly discuss the constraints that we can put on our modelUng. 

• The size R of the emitting region must be compact, in order 

to vary with the observed short timescale. In the adopted model, 
this corresponds to a location of a jet dissipation region Pdiss ~ 
ctvar5/[(l + z)ip] ~ 7 X 10^^tvar,6h('5/20)/V'-i cm, if the jet 
opening angle is = O.l^"-!. 
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Figure 4. The SED of 3C 454.3 at 6 selected epochs: 6 and 27 Nov. 2009, 1, 2, and 3 Dec. 2009 plus the previous "tow" 7-ray state (see text). The different 
SED are labelled. We also show the result of our modelling, including the accretion disk component, its X-ray corona contribution, and the infrared (IR) 
emission from the torus (dashed black lines). For the SED of Nov. 6 we show also, for illustration, the contribution of the EC and SSC components. Note that 
both are necessary to explain the X-ray spectrum, according to our modelling. 



• The short variability timescale and the large disk luminosity 
imply 7?diss < ^Jblr: the dissipation occurs within the BLR. 

• If so, the relevant radiation energy density is pro- 
vided by the BLR photons, with energy density f/eLR ~ 
/BLRl/dr^/(47ri?BLRc)- Setting /blr = 0.1 and using /Jblr = 
10"L^^45 cm, the resulting [/blr ~ rV(127r), i.e. the depen- 
dence on Ld and 7?blr drops out. As a consequence the Lorentz 
factor of the electrons cooling in a time R/c is {Ub and U'^y-a are 
the magnetic and the synchrotron energy densities): 

3moC^ 



(r/20)27?l«[l + [Ub + [/s'yn)/[/^,LR] 

Therefore the radiative cooling is efficient for all but the lowest 
energy electrons. 

• The bulk Lorentz factor F and the Doppler factor 5 are con- 
strained by the observed superluminal motion at the VLBI scale. 
Jorstad et al. (2005) find different components moving with dif- 
ferent apparent velocities /3app, from a few to more than 20, for 
observations performed between 1998 and 2001, resulting in dif- 
ferent bulk Lorentz factors (from 10 to 25) and different orienta- 
tions (viewing angles from 0.2 to 4 degrees), resulting in different 
Doppler factors (<5 from 14 to 30). More recent observations by 
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Lister et al. (2009) found /3app = 14.2 ± 0.8, consistent with the 
previous findings. 

• The soft slope of the 7-ray spectrum and the hard slope of the 
X-ray spectrum constrain the peak of the high energy component 
hfc of the SED to lie below 100 MeV (but close to this value). For 
the EC process, the peak is made by the most relevant electrons 
scattering the Ly~Q seed photons. They must have a Lorentz factor 
7peak given by 

2 _£S 

1 + Z 

>-^~210(^^^J (— j (6) 

This implies that the emission at and above the 7-ray peak is in fast 
cooling (i.e. that the cooling time is shorter than ii/c). 

• Since the bolometric luminosity is dominated by the 7-ray 
emission, this limits the injected power in the form of relativis- 
tic electrons. The fast cooling regime ensures that almost all the 
injected power is converted into radiation. We then have P/^j ~ 

• In general, the value of the magnetic field in the emitting re- 
gion is derived by the level of the synchrotron emission. However 
the peak of the synchrotron spectrum, and hence its bolometric lu- 
minosity, lies in the unobserved mm-far IR frequency range. A 
more robust constrain comes from the soft X-ray band, where we 
have the contribution of the EC and SSC components. Their spec- 
trum is different (the SSC is softer). Increasing the magnetic field 
B increases the SSC component. The B-field is then found by re- 
producing the observed flux and shape of the X-ray spectrum. 

• The observed soft shape of the 7-ray and optical-UV spec- 
tra (once accounting for the contribution of the accretion emission, 
important in the low states) constraints the index S2 of the energy 
distribution of the injected electrons. Since the cooling is almost 
complete (i.e. almost all electrons cool in i?/c), the shape of the 
emitting distribution, below 7pcak, is oc 7^, almost independent 
of si (if si < 2). On the other hand, the details of the curvature 
around 7poak and the distribution below 7cooi do depend on s\. 
This index is then found by reproducing the X-ray spectrum (and 
the 5wi/f/BAT point). 



5 RESULTS 

Fig. |4] shows the results of our modelling for the 6 simultaneous 
SEDs for different epochs, while Tab. |4] lists the used parameters. 
Fig. [5] shows the same SEDs and models, but extending the fre- 
quency axis to the radio band and showing some archival data, to 
put the SEDs analysed here in contest. 

The model successfully reproduces the data from the optical 
to the 7-ray band, using the guidelines explained in i]4.1| The large 
variations of the 7-ray luminosity are obtained by mainly varying 
the power injected in relativistic electrons (by a factor 10 from the 
"/ow" to the highest state, and by a factor 3 in the restricted pe- 
riod from Nov 6 to Dec 2), and the bulk Lorentz factor (from 15 
to 20). These are the main changes. Besides them, there are minor 
changes in the location of the dissipation region (hence, in the size 
of the emitting region) by less than a factor 1.4, and in the parame- 
ters of the injected electron distribution (i.e. 7b changes by a factor 
~ 2, and also si changes somewhat). These changes are required 
to fit the X-ray spectrum, while the changes in the S2 parameter 
are necessary to account for the (small) changes of the 7~ray and 
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Table 4. For all modefs we have a.ssumed a viewing angle d = 1.8° and 
a bolometric lumino.sity of the accretion disk idisk = 6.75 X IC*^ erg 
s~^. The black hole mass is assumed to be M = 5 X 10** Mq, so the 
disk luminosity is the 90% of the Eddington value, and the mass accretion 
rate is Min = 14.78M0 yr~i (assuming an efficiency r\ = 0.08). The 
size of the emitting blob, for the assumed distances Pdiss from the black 
hole, is always i/iPdiss, with ■)/> = 0.1. We assume that the 10% of the disk 
luminosity is reprocessed by the BLR and re-emitted as broad fines at a 
distance J?blr = 8.2 x 10^^ cm for all models. Anotfier 30% of the disk 
luminosity is reprocessed and re-emitted as IR radiation by a torus, located 
at Pin ~ 2 X 10^^ cm. The minimum observed variability timescale ivar 
is in hours and is defined as tvar = -R(l + 2:)/ (c(5). Powers are in units of 
erg s~l. Size and distances in units of 10^^ cm. Magnetic field B in Gauss. 
Outflow mass rates Mout in Mq yr~ ^ . 



the optical-UV spectra. Although minor, these changes have a rel- 
atively strong impact on the total jet power, because 7b and s\ con- 
trol the total amount of electrons present in the source, thus also 
the number of cold protons, if we assume one proton per electron. 
Indeed, the maximum value of Pp is obtained on Dec. 3, not at the 
maximum of the 7-ray flux (occurring at Dec. 2, when there is also 
the maximum value of Pi'nj). 

As explained in i ]4. 11 the value of the magnetic field is mainly 
derived to adequately fit the X-ray spectrum, rather than the syn- 
chrotron optical-UV spectrum and flux level. Fig.|4]shows the con- 
tributions of the SSC and EC components for the SED of Nov 6. 
They intersect at ~1 keV, with the SSC dominating below, and the 
EC above. A similar decomposition has been adopted in Vercellone 
et al. (2010) for another epoch (see their Fig. 19, showing also the 
contribution of the accretion disk radiation to the optical-UV flux). 

Fig.|5]shows that the during the optical flare of 2005 the opti- 
cal flux reached a brighter state than in Nov-Dec 2009. There was 
no 7-ray facility in orbit in 2005, so we do not know the 7-ray flux 
during the optical flare, interpreted by Katarzynski & Ghisellini 
(2007, see also Ghisellini et al. 2007) as a dissipation event occur- 
ring relatively close to the black hole, where the magnetic energy 
density can be of the same order of the radiation energy density as 
seen in the comoving frame. In this case Lsyn ~ i-EC ~ ^7, and 
we do not need an extraordinary increase of the total power budget 
to explain this extraordinary optical flare that can occur without a 
corresponding large increase of the 7-ray flux (and the bolomet- 



10 Bonnoli et al. 



BO 
U 



:^ 

o 
J 



-10 



-11 



-12 - 



1 1 1 1 1 1 

. 3C 454.3 


1 1 1 1 1 1 1 1 1 ^ 

2/12 


" z = 0.859 








2005 flare 


3/ 1 27 / A 




















If 




M M '/I low W 








^ # SAX 




■ 1 
















- s° 1 

^ 1 




I , , , 1 , , , , 





49 



48 ^ 



47 ^ 

O 



46 



10 



15 20 
Log V [Hz] 



25 



Figure 5. Same as Fig.|4] but extending to radio frequencies and sliowing 
also archival data (in light grey), including the optical fluxes achieved dur- 
ing the 2005 optical flare, the 5-6 June 2000 BeppoSAX spectrum (Tavec- 
chio et al. 2002) and the EGRET spectrum of Jan 1992 (Nandikotkur et al. 
2007). Note that i) during the optical flare of 2005 the optical flux reached 
a brighter state than in Nov-Dec 2009; ii) the low Fermi/LAI state is at a 
lower level than detected by EGRET and iii) that the X-ray flux and spec- 
trum detected by 5vv!yf/XRT for the "low" state superpose exactly to the 
BeppoSAX data. 



ric one; in this sense the model can be thought as "economic", i.e. 
assuming the lowest possible power budget). On the contrary, the 
2009 7-ray flare does change the bolometric luminosity by a large 
factor, without inducing a correspondingly dramatic increase of the 
optical flux. The results of the modelling presented here explain 
this behaviour by having a larger power dissipated in regions more 
distant than in the 2005 flare, with relatively smaller magnetic field 
and larger bulk Lorentz factors. 

Fig.HJalso shows that the "low" FermifLAT state is lower than 
the 7-ray flux detected by EGRET in Jan. 1992 and that the X-ray 
flux and spectrum detected by Swift/XKT for the "/ow" state super- 
pose exactly to the BeppoSAX data of 5-6 June 2000 (discussed in 
Tavecchio et al. 2002). 



5.1 7-ray vs X-ray and 7-ray vs optical flux correlations 

In §3 we showed that the 7-ray flux varies almost quadratically 
with the optical and X-ray fluxes. In the framework of synchrotron 
-I- External Compton model, modest variations of the synchrotron 
flux accompanied by large variations of the EC flux imply that the 
magnetic field must vary oppositely with the observed bolometric 
power (i.e. the 7-ray luminosity, in this case). Another way to have 

more than linear variations (but not quadratic) is to vary F, that 

3/2 

would result in Leg oc Lsyn (this relation, valid considering fre- 
quency integrated luminosities, remains valid in restricted energy 
ranges if the EC and synchrotron peaks do not move). 

Fig.[6]illustrates the specific way of our modelling to have the 
observed Lj oc L^pt correlation. The figure shows how B, T and 
/?diss depend on the observed L-,, and how B depends on the prod- 
uct ri?diss- We can see that B inversely correlates with L^, while F 
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Figure 6. Top: the magnetic field B (left) and the bulk Lorentz factor F 
(right) as a function of the 7-ray luminosity Lj for the 6 different SED 
plotted in Fig.|4] Bottom left: the magnetic field B as a function of FiJ^igg . 
Bottom right: -Rdiss a function of Lj . 



and i?diss, instead, correlate positively. These makes B to be almost 
perfectly proportional to (FTidiss)"^- Since Pb oc [PdissFB]^, we 
have that the Poynting flux at Rdiss is approximately constant for 
all the considered states. 

To summarise: the large variation in the 7-ray flux accompa- 
nied by much more modest variations of the optical are explained 
by a magnetic field that decreases when the observed increases. 
This is accomplished, in our model, by having larger i?diss and 
F for larger and jet powers, and by having instead a quasi- 
constant Poynting flux in the emission region. By the same argu- 
ment we can explain the F^-Fx correlation, since the contribution 
of the SSC component to the X-rays becomes relatively more im- 
portant for lower states, making the SSC X-ray flux to decrease 
much less than the 7-ray one. 

However, the quasi constancy of Pb may be partly a coinci- 
dence, and not the outcome of a fundamental process. This is be- 
cause it is very likely that the Poynting flux Pb,o at the start of 
the jet is the prime cause of the acceleration of the jet itself. In 
this case, since we find that the total jet power does change, Pb,o 
should change as well, and yet have the same value when arriving 
at Pdiss- What can happen is that a stronger Pb.o implies a larger 
final F, achieved at a larger distance z (i.e. we can have F oc 2^''^ 



as in Vlahakis et al. 2004). This means that a larger Pj, 



jet 



B,0 IS 



achieved, but since both Pdiss and F are larger, Pb at the dissipa- 
tion region may vary much less than Pb,o ~ Pjct oc L-,. 



5.2 The power of the jet 

It is interesting to investigate if the extraordinary 7-ray flare cor- 
responds to a real enhancement of the power carried by the jet or 
if originates by an enhanced efficiency in converting the jet power 
into radiation, or if is due to an increased Doppler factor 5, making 
the observed flux to increase even if the jet power remains constant. 
The nearly model-independent quantity connected to the jet power 
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Figure 7. The observed 7-ray luminosity (top panel) the power 
spent by the jet to produce the radiation we see (mid panel) and the total 
jet power Pjet (botton panel) as a function of the accretion disk luminosity 
Ljj. The filled circles correspond to the 6 states of 3C 454.3 analysed in this 
paper, while the empty squares con'espond to all the FSRQs detected in the 
first 3-months of the Fermi/LPH all sky survey analysed by Ghisellini et al. 
2010. 



is Pr, the power that the jet spends to produce the radiation we see. 
Using U^.^^ = L' /{A'kB?c), we can re-write Eq.|4]as 



L 



y 



L 



(7) 



where L is the total observed non-thermal luminosity (i' is in the 
comoving frame) and Ul.^^ is the radiation energy density produced 
by the jet (i.e. excluding the external components). The last equality 
assumes 6^ ~ l/F. This quantity is almost model-independent, 
since it depends only on the adopted 5, that can be estimated also 
by other means, namely superluminal motions. In Tab.|4]we list the 
values of Pr for the different states. The maximum is on Dec 2, the 
day of the maximum L-y, when Pr = 5.2 x 10*^ erg s~^. This value 
can be compared with the accretion luminosity, that we assumed 
constant during the entire period and equal to Ld ~ 6.75 x 10'*'' 
erg s~^. The power Pr can be considered as a very robust lower 
limit to the total jet power Pjot = Pb + Po + Pf 



Pjet is used to produce the radiation we see (i.e. Pj 



If the entire 
jet = Pr), then 

the jet would necessarily decelerate and there would be no VLBI 
superluminal knots, nor any extended radio emission. 

Tab.|4]lists also the value if Pb and Po. They are smaller than 
Pr. The fact that Pc < Pr may seem strange at first sight, since the 
radiation is produced by the electrons, so how can they have less 
power that what they transform in radiation? The answer is that the 
radiative cooling time is shorter than R/c, so electrons need to be 
injected continuously, at least for a time R/c. To estimate Pc we 
count electrons (and their energy) forming the N{'y) [cm~^] en- 
ergy distribution after solving the continuity equation, not the ones 
injected. An example can clarify this point: assume to inject the 
same energy distribution (3(7) [cm^"^ s^*]. In case A the radiative 
cooling 7 is strong, in case B is weak. The particle density N{-^) is 
always proportional to [J Q(7)d(7)]/7, therefore in case A A''(7) 
is smaller than in case B. This despite the fact that in case A we 
have produced more radiation, since the cooling is stronger. 

This implies that there must be a "reservoir" of power able to 
energise electrons that in turn produce Pr. The simplest solution is 
to assume that this "reservoir" is provided by protons. The listed 
values of Pp assume one proton per electron. The implication is 
that electron-positron pairs cannot be energetically dominant, even 
if we cannot exclude that they outnumber primary electrons (but not 
by a large factor). For a more detailed discussion about this point 
see Sikora & Madejski (2000) and Celotti & Ghisellini (2008). We 
summarize here the arguments in Celotti & Ghisellini (2007). If 
pairs are dynamically important, they must outnumber protons and 
also the emitting leptons by a large factor. Although possible in 
principle, this solution poses the problem of the origin of this large 
number of pairs. If they have been created in the same emitting 
region (e.g. by 7-7 collisions), then we should see i) a cutoff in 
the spectrum, and ii) some reprocessed radiation in X-rays (since 
the pairs are born relativistic). The level of this reprocessed radia- 
tion must be of the same order of the absorbed luminosity (i.e. it is 
large and should be well visible). If, instead, the pairs are created 
from the start, we can calculate how many we need. In powerful 
sources (as 3C 454.3) this number corresponds to an initial pair op- 
tical depth greater than unity. If initially cold, the pairs quickly an- 
nihilate, and the surviving ones are not enough (to be dynamically 
iinportant). If they are hot, they emit. At the start of the jet, close 
to the accretion disk, the radiation energy density and the magnetic 
field make them to cool very rapidly. So they become cold, and 
annihilate (besides producing radiation we do not see). 

Assuming then that Pp is a good proxy for the real Pjot we 
see that during Nov and Dec 2009 the jet power varied only by a 
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factor of 2, while varied by a factor ~10. This implies that the 
jet became more efficient in transforming its power into radiation, 
rather than becoming more powerful. For the "/ow" state, instead, 
the estimated value of Pp ~ 10^^ erg s~^ is a factor 10 less than 
in Nov 27, with Lj being a factor 4 weaker. To summarise: the 
total excursion (from "/ow" state to Dec 3, 2009) of Pp is factor 20, 
not very different to the total amplitude of L-y (factor ~30). But 
restricting the period from Nov 6 to Dec 3 the jet power is quasi- 
constant, while L-y varied by a factor 10. 



5.3 Jet power vs disk luminosity 

The disk luminosity is Ld = 6.7 x 10*® erg s~^, comparable to Pr. 
The jet power must be greater than Pr and therefore greater than 
Ld. This is one of the most important outcomes of having followed 
3C 454.3 during its major 7-ray flare. 

Tab. |4] lists the value of the outflowing mass rate, estimated 
through Pp — VMoutf? ■ Since V varies moderately (between 15 
and 20) while Pp varies by a factor 10, we have that also ALout 
varies by an order of magnitude considering the entire time span, 
and only by a factor 2 in Nov - Dec 2009. In this period it is of the 
order of 1 solar mass per year. We can compare this value with the 
accretion mass rate Min, that we can derive through Ld — rfMinC? ■ 
Assuming 77 — 0.08, a disk luminosity Ld — 6.75 x 10*® erg s~^ 
gives Afin ~ I4.8AL0 yr~^, about a factor 10 greater than ALout. 

If the disk emits at a constant level, with a constant ALin, we 
are forced to conclude that the link between the accretion rate and 
the jet power is not determined by ALin, or, rather, that this is not 
the only important quantity in producing Pjet. For the ensemble 
of bright FSRQs detected by Fermi, in fact, the disk luminosity 
does correlate with Pjct, even taking out the effect of a common 
redshift dependence in the two quantities (see Fig. |7] and the dis- 
cussion in Ghisellini et al. 2010). Since the spin of the black hole 
of 3C 454.3 is constant, the likely quantity that modulates Pjct in 
Nov-Dec 2009 is the value of the magnetic field in the vicinity 
of the black hole horizon. If we assume that Pb,o = Pjct at the 
Schwarzschild radius Ps, the magnetic field must be 



Bo = B(Rs 



ir,5 pl/2 
-^jct,48 



5 X lO^ALe 
M 



1/2 



G 



(8) 



In turn, we can check if the magnetic energy density Bq/{8it) can 
be equated to pc^, the mass energy density of the accreting matter 
close to the Schwarzschild radius. We can write 



pc 



9 H 

A-kRI — pV^ ^ 
its 

Afinc" 

A-kRI{H/Rs)Vb. 



(9) 



47rr;Pi(7L/Ps)1^R 

pc^ _ _Ld 1 

Pg/(87r) " Pjct MH/Rs)^R 

where Vr is the radial infalling velocity and 2H is the total height 
of the disk at the Schwarzschild radius. The above equation shows 
that if pc^ ~ Bo/(87r) and Pjct/Ld ~ 10, then Pr{H/R) ~ 0.3, 
i.e. not an unreasonable value. 



5.4 Is 3C 454.3 exceptional? 

We finally ask if 3C 454.3 is exceptional, or instead if there are 
other blazars reaching comparable values of and of L-^/Ld. 



To answer, we can compare the flaring state of 3C 454.3 with all 
blazars of known redshift of the first 3-months FermifLAT all sky 
survey (Abdo et al. 2009c), analysed in Ghisellini et al. (2010), as 
done in the top panel of Fig.|7] We alert the reader that in Ghisellini 
et al. (2010) we have considered the average 7-ray flux. Had we 
considered the peak flux values there would be a few FSRQs of 
L-y even greater than 3C 454.3 in high state (one example is PKS 
1502+102, Abdo et al. 2010a). In this respect we can conclude that 
3C 454.3 is not the most extreme blazar: its unprecedented 7-ray 
flux is due to its relative vicinity in comparison to other FSRQs. 
The same occurs in the planes Pr-Ld and Pjct-Ld (mid and bottom 
panels of Fig.|7](: there are other sources whose 3-months averaged 
Pr and Pjct are comparable to 3C 454.3. 



6 CONCLUSIONS 

The study of the strong 7-ray flare of 3C 454.3, together with the 
good coverage at X-ray and optical-UV frequencies provided by 
Swift allowed to investigate several issues about the physics of the 
jet of this blazar. In Tavecchio et al. (2010) we analysed the be- 
haviour of the 1.5 year FermifLAT light curve, finding episodes 
of very rapid variations, with time scales between 3 and 6 hours 
both during the rising and the decaying phases. This implies a very 
compact emitting region, suggesting that the dissipation zone is not 
too far from the black hole, and that cooling times are shorter than 
the variability time scales. In turn this suggests that the dissipation 
region lies within the BLR, at about one thousand Schwarzschild 
radii. 

We have revisited the estimate of the mass of the black hole 
of 3C 454.3, finding a smaller value than found by Gu, Cao and 
Jiang (2001). Since in our model the mass of the black hole sets the 
scales of the disk+jet system, a smaller mass helps to have a more 
compact emitting region, that can vary on shorter time scales. 

The optical. X-ray and 7-ray fluxes correlate. This supports 
one-zone models. The 7-ray flux varies quadratically (or even 
more) with the optical and X-ray fluxes. By modelling the opti- 
cal to 7-ray SED with a one-zone synchrotron+inverse Compton 
leptonic model we can explain this behaviour if the magnetic field 
is slightly fainter when the overall jet luminosity is stronger. 

The power that the jet spent to produce the peak 7-ray lumi- 
nosity is of the same order than the accretion disk luminosity. Al- 
though the jet power correlates with the accretion luminosity con- 
sidering the ensemble of bright FSRQs detected by Fermi, 3C 454.3 
probably varied its jet power while maintaining a constant accretion 
luminosity. This implies that the modulation of the jet power may 
not be due to variations of the accretion rate, but is probably due to 
variations of the magnetic field close to its black hole horizon. 
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